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Abstract. Fluxon transmission through several impurities of different strength and type (i.e., microshorts 
and microresistors) , placed in a long Josephson junction is investigated. Threshold pinning current on the 
impurities is computed as a function of the distance between them, their amplitudes and the dissipation 
parameter. It is shown that in the case of consequently placed microshorts or microresistors, the threshold 
pinning current exhibits a clear minimum as a function of the distance between the impurities. In the case 
of a microresistor, followed by a microshort, an opposite phenomenon is observed, namely the threshold 
pinning current exhibits maximum as a function of the distance between the impurities. 

PACS. 03.75.Lm Josephson vortices - 05.45.Yv Solitons - 74.50.+r Josephson effect 
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1 Introduction and the background 

The dynamics of magnetic flux propagation in a long Joseph- 
son junction (LJJ) is a subject of increasing theoretical 
and practical interest [HE]- Magnetic flux quantum in a 
LJJ is a soliton (also known as fluxon) governed by the 
well-known sine-Gordon (SG) equation. A convenient way 
to prepare a junction with the required properties is to in- 
stall various inhomogcncitics into it. Up to now substantial 
work has been devoted to the study of the fluxon motion 
in the LJJs with point-like impurities. The interaction of a 
fluxon with a single impurity became a textbook example 

On the other hand, the phenomenon of resonant tun- 
neling of an electron through a double-well structure is 
well-known in quantum mechanics [4] . A natural question 
arises: what is an analog of the quantum-mechanical res- 
onant tunneling in the fluxon dynamics? Resonant soli- 
ton transmission has been investigated in detail for non- 
dissipative systems [5J|S] and complex resonant behaviour 
has been reported. However, fluxon dynamics in a LJJ 
cannot be considered without taking into account dissipa- 
tive effects, which are a consequence of the normal elec- 
tron tunneling across the insulating barrier. As a result, 
transmission in a LJJ with constant bias and dissipation 
can yield only two scenarios: fluxon transmission or fluxon 
pinning on the impurities. And, consequently, the trans- 
mission ratio can attain only two values: zero or unity. 
Therefore the attention has to be turned toward other 
characteristic quantities, especially the minimal bias, nec- 
essary for the fluxon pinning on impurities. 

The present paper aims to investigate fluxon transmis- 
sion through several (two or more) point-like impurities: 
microshorts, microresistors or a combination of both. Of 
particular interest is dependence of the threshold pinning 



current on the distance between the impurities and their 
amplitudes. 

The paper is organized as follows. In Section [2] we 
present the model and the basic equations of motion. In 
the next section we describe the methods of the analysis 
of the equations and motion and study the fluxon trans- 
mission through two microshorts, two microresistors and 
a microshort and a microresistor as a function of their 
amplitudes and distance between them. Discussion of the 
obtained results and final remarks are given in Sec. [U 



2 The model 

We consider the long Josephson junction (LJJ) subjected 
to the external time-independent bias. The main dynami- 
cal variable is the difference between the phases 02{x,t) — 
#i(x, t) = 4>(x, t) of the macroscopic wave functions of the 
the superconducting layers of the junction. The time evo- 
lution of the phase difference is governed by the perturbed 
sine-Gordon (SG) equation: 

4>tt~ 4>xx + sin 4> = ef [(f), <j> t ]x] , 

m 

ef[(f>,(j>t;x} = -acpt - 7 - sin</>^ /J, n 8(x - a m ). (1) 

n=l 

In this dimensionless equation spacial variable x is nor- 
malized to the Josephson penetration depth Aj, the tem- 
poral variable t is normalized to the inverse Josephson 
plasma frequency wj 1 !?]. Here the bias current 7 is nor- 
malized to the critical Josephson current of the junction 
and the dimensionless parameter a describes dissipation. 
It is supposed that there are N impurities in this junc- 
tion, positioned at the points x = a„, n — 1, 2, ...,N, 
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a\ = < a 2 < ... < a m , with \x n being "strength" or am- 
plitude of the nth impurity. The impurity is a microshort 
if \i n > and a microresistor if /i„ < 0. 



3 Fluxon transmission 

A standard tool for analyzing the fluxon dynamics in Joseph- 
son junctions is the McLaughlin-Scott perturbation theory 
|3J. Also, direct numerical integration^ of the perturbed 
SG equation (pj) will be performed to check the validity of 
the analytical approximation. We are going to solve the 
problem for the idealized case of an infinite junction with 
free ends boundary conditions, however, in actual simu- 
lation a sample with length that significantly exceeds the 
fluxon size will be used. 



3.1 Perturbation theory and collective coordinates 

Using the perturbation theory, one obtains in the first 
order the evolution equations for the fluxon parameters, 
i.e., its center of mass X and fluxon velocity v: 

v = ^-(l-v 2 ) 3 / 2 -av(l-v 2 ) + 



In the case of one microresistor fi < 0, the threshold bias 
has been defined in [8 as 



1 - v 2 



N 



^2n n g[X -a(n- l),v] , 



n=l 
V 



N 



X = v - - Vn[X - a(n - l)]g[X - a(n - 1), «],(3) 



g(x,v) 



tanh (X/Vl-v 2 ) 
cosh 2 (X/VT^ 7 ^ 2 ) 



For the sake of simplicity in the following only equidistant 
impurities will be considered, i.e., a n = a, n = 1, 2, . . . , N. 
Also, only positive values of bias 7 will be considered. The 
case of one impurity (N = 1, fi\ = fi) has been discussed 
in detail in [3JH]. There exist two characteristic values of 
the bias current, j c = 4-^/3^/(9^) and 7*^,-, 7 C > Jthr- If 
7 > 7c, the pinning on the impurity is not possible and 
only one attractor that corresponds to fluxon propagation 
does exist. In the interval jthr < 7 < 7c two attractors 
exist: one corresponds to fluxon pinning on the microshort 
and another one to fluxon propagation. If 7 < jthr, the 
only possible regime is fluxon pinning on the impurity. It 
has been shown [3j[8] that there exists a threshold value 
of the dc bias, which can be approximated as 



Ithr 



a 



7T 



fi 2 [1 - 2aln2] 



(4) 



1 In the numerical simulations, the space will be discretized 
as x — » nh, so that the continuous variable 4>(x, i) ~ <f>(nh, t) 
becomes a discrete set of variables 4>-n{t), and the second space 
derivative becomes <f) xx (x,i) ~ [(j> n +i (t) -2<f> n (t) + (f> n - 1 (t)]/h 2 . 
The resulting set of the second order ODEs on (f> n (t) will be 
solved using the 4th order Runge-Kutta scheme. The delta 
function is approximated as S(x) ~ 8„,o/h where S m ,n is Kro- 
necker's S symbol. 
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-am(a^/2/M) 



(5) 



In the non-relativistic limit (v <§C 1) the system ([2])-(j3]) 
can be rewritten as a Newtonian second order ODE for 
the particle of mass 8 (see Refs. [H[8]): 

dU 

&X + 8aX = -—, (6) 

U(X) = -2tt 7 A + U (X) = -2n-fX + (7) 
m 



+2 



^ cosh 2 [X 



1)] 



7 > 



In the case of strongly separated impurities (n> 1), the 
potential U(X) has 2N extrema [they approximately co- 
incide with the fixed points X = X 2 k, k = 1, 2, . . . , N, of 
equations ([2])- (J3j)] where each pair (a minimum and a max- 
imum) is associated with a certain impurity. If there is an 
impurity at X = a(k— 1), the minimum at X = X^h-i al- 
ways comes before the maximum at X = X^k- Microshorts 
are repelling impurities, thus the fluxon that arrives from 
X = — 00 decelerates when approaching it and accelerates 
after passing the impurity until the fluxon velocity reaches 
the equilibrium value 



, 4a 
1+ I — 

7T7 



-1/2 



(8) 



Microresitors are attractive impurities, and, as a result, 
the fluxon accelerates before approaching the microresis- 
tor and slows down to the equilibrium velocity {5J after 
being released from it. Decrease of the distances between 
the impurities a causes disappearance of some of the ex- 
trema via inverse pitchfork bifurcations. The systematic 
phase plane analysis of equations ©-fS]) for the case of two 
microshorts has been performed in [9] for the SG equation 
and in [10] for the double SG equation. In those papers 
the behaviour of the fixed points of the system ([2|)-(l3]) has 
been studied as a function of the distance between them, 
a. 

Our aim is to determine the threshold current jthr = 
1thr(a, {Hn}\ ] N) as a function of the distance between 
impurities and their amplitudes. In the case of one impu- 
rity, 7t/ l r(A t i 1) obviously does not depend on the distance 
a. It is described approximately by equations (J4| and (O 
for the microshort and microresistor, respectively. Some 
general statement can be made before one proceeds to 
specific cases. Two important limits should be mentioned. 
One case corresponds to impurities being separated by 
the distance much larger the fluxon size. Then the trans- 
mission will be governed by the fluxon interaction with 
each individual impurity. In the opposite limit (a — > 0) 
the power of all impurities adds up. The effect of the both 
limits on the threshold current can be written as follows 



lthr{a,{Hn}i ;N) = 



0. 



(9) 
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In the subsections below the transmission through im- 
purities of different polarities (e.g, ji n < and fj, n > 0) 
will be considered. It appears that for N > 2 the analyti- 
cal treatment of equations ©-([3]) is virtually not possible 
even in the non-relativistic case, especially when none of 
the limits, described by the equation ([9]), hold. Therefore 
equations Q-© are going to be solved numerically. 

3.2 Transmission through two microshorts 

Consider first the case of two microshorts (N — 2, \x\p. > 
0). The problem is tackled in the following way. The fluxon 
approaches the system of two microshorts from X — — oo 
with the equilibrium velocity Voo, given by equation ([8]). 
Evolution of the system ©-([HI on the phase plane (X,v) 
is shown in Figure [TJ Depending on the strength of the 
bias, three scenarios are possible: trapping on the first mi- 
croshort (curve 1); trapping on the second microshort, if 
the external bias is a bit larger (curve 2); or transmission 
(curve 3). If the microshorts are too close to each other, 




Fig. 1. Phase space trajectories of equation ([2} for N = 2 
microshorts with /ii = 0.4, fi2 = 0.6 and distance a — 3 be- 
tween them. Curve 1 corresponds to 7 = 0.04, curve 2 - to 
7 = 0.05 and curve 3 - to 7 = 0.06. Blue curve shows results 
of direct integration of the equation {TJ at 7 = 0.06 the using 
Runge-Kutta method. The fluxon center of mass is defined as 
X(t) = f^^x&xdx and its velocity as v(t) = dX(t)/dt. Dis- 
sipation for all cases is a = 0.1. The inset shows the effective 
potential U(X) [equation 0] for /xi = 0.4, fi 2 = 0.6, 7 = 0.05, 
a — 1 (curve 1), a — 1.7 (curve 2), a = 2.5 (curve 3) and a — 4 
(curve 4). Red vertical bars denote locations of the microshorts 
and the lengths of the bars is proportional to the microshort 
amplitudes, [lip. 

trapping on the second microshort does not happen (see 
references [91H0] for details). We note that direct numeri- 
cal simulations of the perturbed SG equation ([TJ (curve 4 
of Figure[TJ are in good correspondence with the trajecto- 
ries of the system ©. The oscillations after the collision 



with the microshort can be attributed to the fluxon radia- 
tion (not accounted by the first order perturbation theory) 
and errors in determination of the fluxon center. 

The systematic evaluation of the threshold current "f t hr 
as a function of the distance a for different values of fi± 
and H2 is shown in Figure [5J The resonant nature of the 
dependence of jthr on a for < 1x2 can be observed 
clearly. While in the respective limiting cases it satisfies 
equation ©, a resonant value a = a r appears, at which 
the threshold current attains its minimal value. 



Y 




2 4 6 8 10 

a 



Fig. 2. Threshold pinning current as a function of the distance 
between the microshorts for /ii = 0.6, fj,2 = 0.4 (A); fii = 
fi 2 = 0.5 (o); m = 0.4, H2 = 0.6 (+); yui = 0.3, fi 2 = 0.7 (x); 
/ii = 0.2, H2 — 0.8 (V). Results, obtained via direct numerical 
integration of the equations of motion for fii — 0.4, [ii = 0.6 
are shown by 0. 



The explanation of the resonant transmission can be 
done on the following qualitative argument. The analysis 
of the phase portraits in Figure [TJ shows that after be- 
ing released from the microshort, the fluxon accelerates in 
order to regain its equilibrium velocity Voq . This accelera- 
tion occurs in such a way that for some short interval the 
fluxon velocity exceeds the equilibrium value Voq. There- 
fore the fluxon has kinetic energy, which is larger than 
it was while approaching the microshort from X — — 00, 
and consequently it has enough energy to pass the mi- 
croshort with the amplitude larger than Obviously, 
the best transmission would take place if a slightly ex- 
ceeds |X2 — X\\. The estimation of the resonant distance 
a r can be made from the analysis of the fluxon dynam- 
ics in the non-relativistic limit, given by equations (|rJj)- 
dTj) . According to these equations the fluxon can be com- 
pared to the particle that slides down along the poten- 
tial U(X) = U(X -> ±00) 2ir-yX. Depending on the 

value of 7, it can be trapped in one of the wells of this 
potential (shown in the inset of Figure [JJ . If the distance 
between microshorts is small enough, it can be consid- 
ered as one microshort with the renormalized strength 
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p(a) = fi% + cosh 2 a. The trapping can occur at the 
only existing minimum X\ (see curve 1 in the inset of 
Figure [I]) as shown by the trajectory 1 of Figure [TJ As 
a increases, the potential barriers separate and a local 
minimum X3 appears, as shown by curves 2 — 4 in the 
inset of Figure [I] If a new minimum appears, the trap- 
ping can occur also at the second microshort, as shown 
by the trajectory 2 of Figure [TJ Plots of the potential 
U (X) clearly demonstrate that the shape of the barrier 
will be optimal when the minimum at X — X3 is quite 
shallow. Since the half- width of the function cosh _2 (X) 
is of the order of unity, it is expected that optimal sepa- 
ration of barriers occurs at a ~ 2. Numerical evaluation 
of a r confirms this estimate: a r — 1.94 (for /ii = 0.4, 
/j 2 = 0.6); a r = 2.24 = 0.3, /i 2 = 0.7) and a r = 2.36 
(for iii = 0.2, /12 = 0.8). If /ii > /12, the transmission 
scenario is always determined by the first microshort and 
the trapping occurs only at X = X\. Therefore the de- 
pendence jthr on a is monotonically decreasing as shown 
in Figure [2] for fii = 0.6, [i% = 0.4. 

It would be of interest to compare how the threshold 
pinning current depends on the dissipation parameter a 
and the ratio of the microshort amplitudes /ii and /i2- 
Since the resonant distance a r weakly depends on fii^ and 
the pinning current depends strongly on the dissipation 
constant, it is convenient to normalize ^thrifl, fii, fj,2', 2) to 
the pinning current on the strongest microshort 
max [ythrifJ'V, 1), 7thr(^2] 1)]- In Figure [3] the dependence 
of the enhancement factor 



77(0,^1,^2; 2) 



7t/ir(a,A t i ) M2;2) 



:[7tfir(Mi;l),7tfcr(/*2;l)] 



(10) 



on the ratio + ^2) for different values of dissipa- 

tion is shown. The value of the distance between the mi- 
croshorts has been fixed to a = 2. Increase of dissipation 
does not change much the resonant values of /j.i.2- How- 
ever, the value of the enhancement factor at the minimum 
decreases significantly. In the inset of Figure [3] comparison 
of the fluxon slowing down on the microshorts is shown 
for different values of dissipation and dc bias. Note that 
the ratio a/7 was kept constant in order to fix the equi- 
librium velocity Wqq. For stronger dissipation the fluxon 
slows down to smaller velocities (compare the black and 
blue curves that correspond to a — 0.1 and a = 0.3, re- 
spectively). Therefore after release from the microshort 
the fluxon can accelerate to greater values of velocity. As 
a result, it has more kinetic energy to pass the second mi- 
croshort. In other words, for larger dissipation one needs 
larger bias, 7. Therefore, the tilt of the potential U(X) 
increases and it smears out the inhomogeneities, created 
by the impurities. 

It should be emphasized that the validity of the per- 
turbation theory approach has been confirmed by the di- 
rect numerical integration of the original perturbed SG 
equation {TJ). In Figure [21 ^ t hr has been computed via in- 
tegration of equation (TTJ) for /ii = 0.4 and [12 — 0.6. It 
is evident that the perturbation theory gives qualitatively 
the same result and the quantitative difference is not very 
large. Similarly, in Figure[3]the results of the numerical in- 
tegration of equation (TTJ with a = 0.3 are given alongside 
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Fig. 3. Dependence of the enhancement factor (|10[) on the 
ratio + /12) a — 0.1 (o), a — 0.2 (o) and a = 0.3 (+). 

Solid lines are used as a guide for an eye. Results, obtained via 
direct numerical integration of the SG equation with the 4th 
order Runge-Kutta method for a — 0.3 are shown by ©. The 
inset shows trajectories in the phase plane (X,v) for a — 0.1, 
7 = 0.06 (black line); a = 0.2, 7 = 0.12 (red line) and a = 0.3, 
7 = 0.18 (blue line). Other parameters are /ii = 0.4, /i2 = 0.6, 
a = 2. 



with the perturbation theory results. A good qualitative 
and quantitative correspondence between these two types 
of results is clearly demonstrated. Therefore the usage of 
the approximation (HI)-© is justified. 



3.3 Transmission through N > 2 microshorts 

Now we extend the results of the previous subsection on 
the case of more then N = 2 microshorts. In Figure [4] the 
dependence of j t hr on a for N = 2,3,4,5 is presented. 
It clearly demonstrates that addition of an extra impu- 
rity to the left from the weakest one decreases further the 
minimum of the threshold current. 

The explanation can be easily seen with the help of 
the effective potential U(X) [see equation 0]. Its shape 
changes significantly when extra microshorts are added. 
Comparing curves 1 and 2 in the inset of Figured] one can 
see that the energy barrier, which the fluxon should cross, 
lowers. Adding yet another microshort further lowers the 
barrier (see curves 3 and 4), so that in the interval < 
X < (N — l)a the potential barrier almost turns into the 
decaying slope which is less steep then — 2-7T7A. Decrease 
of a leads to the gradual raising of this slope (compare 
curves 4 and 5) and consequently to increase of jthr- 



3.4 Transmission through microresistors 

A microresistor is an attracting impurity, therefore the 
fluxon accelerates when approaching it and decelerates 
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a 

Fig. 4. Threshold pinning current as a function of the dis- 
tance between the microshorts for a = 0.1, N = 2, /Lti = 0.4, 
fi 2 = 0.6 (+); N = 3, m = 0.27, (jl 2 = 0.4, ^ 3 = 0.6 (V); 
N = 4, /zi = 0.18, ^2 = 0.27, ^ 3 = 0.4, fi 4 = 0.6 (o) and 
N = 5, fxi = 0.12, ^ 2 = 0.18, ^ 3 = 0.27, ^ 4 = 0.4, ^ 5 = 0.6 
(x). Solid lines are used as an guide for an eye. The inset de- 
picts the effective potential U(X) [see equation 0] for the 
configurations, described in the main figure at 7 = 0.05 and 
for the same values of {/-i n }'- N = 2, a = 2 (curve 1); N = 3, 
a = 1.7 (curve 2); AT = 4, a = 1.6 (curve 3); N = 5, a = 1.6 
(curve 4) and a — 1 (curve 5). 



back to v = Voo after passing through or remains trapped 
if its velocity (and consequently the external bias current) 
is less than the threshold value. The effective potential 
U (X) for a microresistor corresponds to the potential well. 
If two different microresistors are added consequently, the 
fluxon can be trapped on the first or on the second one, 
or, if the bias is large enough, pass through. In Figure 
the phase portraits for the system with N = 2 microre- 
sistors is shown. The change of the shape of U(X) for the 
different distances between the microresistors is shown in 
the inset of Figure [5] 

The computation of the threshold current jthr shows 
that resonant fluxon transmission is possible if fix < fi2 
and does not happen if [i\ > /i 2 (see Figure [6^)- Ex- 
planation of this phenomenon is similar to the case of 
two microshorts. If the microresistors are located very 
close to each other, then their amplitudes add up and 
the fluxon interacts with the microresistor of the ampli- 
tude /i ~ H\ + f/,2- When the impurities start to separate, 
the effective energy barrier which the fluxon should sur- 
mount, lowers (compare curve 1 with curves 2 and 3 in the 
inset Figure [6]). The distance between the wells becomes 
optimal for the best fluxon transmission before they are 
completely separated (compare curves 3 and 4). 

In contrast to the transmission through two microshorts, 
the resonant value a = a r depends strongly on the ampli- 
tudes fii.2- Indeed, for [i\ — —0.7, n% = —0.3 one obtains 
a r — 3.62 and for fi± = —0.9, ji2 = —0.1 the resonant dis- 
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Fig. 5. Phase space trajectories of equations ([111-© for N — 2 
microresistors with /ii = —0.7, /12 = —0.3 and distance a — 2.5 
between them. Curve 1 corresponds to 7 = 0.015, curve 2 - 
to 7 = 0.03 and curve 3 - to 7 = 0.045. Dissipation for all 
cases is a = 0.1. The inset shows the effective potential U(X) 
[equation Q] for — —0.7, /12 = —0.3, 7 = 0.03, a = 1 (curve 
1), a = 2.5 (curve 2), a — 3.5 (curve 3) and a = 5 (curve 4). 
Red vertical bars denote locations of the microresistors and the 
lengths of the bars is proportional to their amplitudes, fJ,i,2- 

tance equals a r ~ 2.35. Comparing curves 2— 4 in the inset 
of Figure [5l one can notice that the fluxon needs enough 
kinetic energy to overcome the second maximum, located 
at X = X4. Obviously, if Xi and X^ are not enough sepa- 
rated, the fluxon will have no time to accelerate in order to 
avoid trapping on the second microresistor. Therefore the 
case of curve 3 is the most optimal one: the height of the 
barrier at Xi is not too large, as compared to the curve 
4 and the distance between X2 and X3 is enough to gain 
velocity, sufficient for the successful passage over the sec- 
ond barrier. These considerations, of course, correspond 
to the situation, when the impurities are not strongly sep- 
arated. Otherwise only the interaction with the first one 
would matter. For the same reason the position of the min- 
imal threshold current, a r , (see Figure Eb) increases with 
decrease of the damping parameter a. Depth of the min- 
imum decays with decrease of a similarly to the case of 
two microshorts because with the stronger bias the fluxon 
can pass through the impurities much easier. 

Putting additional microresistors after /Z2, > (In > 
• • • > > H\ further lowers the critical pinning current 
similarly to the case of N > 2 microshorts, described in 
the previous Subsection. 

3.5 Transmission through a microshort and a 
microresistor 

Finally, we consider the case when two impurities of differ- 
ent polarity (a microshort and a microresistor) are placed 
one after another. If the microresistor is located before 
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Fig. 6. Panel (a). Threshold pinning current as a function of 
the distance between two microresistors for fii — —0.3, /i2 = 
-0.7 {A); m = fi 2 = -0.5 (o); /ij = -0.6, /z 2 = -0.4 (+); /u = 
-0.7, H2 = -0.3 (x); m = -0.8, m = -0.2 (V); jUi = -0.9, 
\i2 = -0.1 (o) and /ti = -0.97, /t a = -0.03 (0). Damping 
coefficient equals a = 0.1. Panel (b). Enhancement factor r\ 
[see equation (|10[) ] as a function of a for /ti = —0.8, fi2 = —0.2 
and different values of damping a: a — 0.3 (+), a — 0.2 (x), 
a = 0.1 (o) and a = 0.05 (V). 



the microshort (/ij < 0, ^2 > 0) resonant enhancement of 
the threshold pinning current does not happen. In Figure 
[7] (panel a) the phase portraits for this case are shown. 
The microresistor is an attracting impurity after which 
the fluxon slows down. On contrary, the microshort is a 
repelling impurity and the fluxon slows down when ap- 
proaching it. Therefore it is obvious that by placing impu- 
rities in such a way one increases jthr as compared to the 
case of each individual impurity. The analysis of the effec- 
tive potential U(X), shown in Figure [7] (panel c), further 
confirms above considerations. The height of the effective 
barrier, which the fluxon should overcome, can be greater 
than the height of the individual barriers, created by the 
individual impurities. If the impurities are very close their 
influences cancel each other and the fluxon interacts with 
the impurity of the strength — |/ti| + (12- The dependence 
of the threshold pinning current ■ythr on the distance be- 
tween the impurities is given in Figured] (panel a). For 
Hi = —fj,2 — —0.5 the microshort and microresistor cancel 
each other for a — 0, therefore the dependence starts at 
zero and increases until it reaches the maximal value and 
then decreases, tending monotonically to the value of "f t hr 
that corresponds to one microshort with fi — 0.5. The de- 
pendence of jthr on a shows an "antiresonant" behaviour 
because it has a maximum at some certain value of a. 
Analysis of the shape ofU(X) from the Figure Eh predicts 
that the worst transmission would occur when the poten- 
tial well and the barrier, created by the microresistor and 
the microshort, respectively, separate from each other far 
enough to create the highest total barrier (see the curve 
3 of Figure |7t) , but not too far (as for the curve 5 of the 
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Fig. 7. Phase portraits of the fluxon dynamics and effective 
potentials for the fluxon dynamics in the case of one microshort 
and one microresistor. Panel (a) shows fluxon trajectories on 
the phase plane (X, v) for a = 0.1, fii = —0.4, fi2 = 0.6, a = 4; 
7 = 0.05 (curve 1), 7 = 0.06 (curve 2) and 7 = 0.07 (curve 3). 
Panel (b) shows fluxon trajectories on the phase plane (X, v) 
for a — 0.1, /ii = 0.2, fi2 = —0.8, a — 3; 7 = 0.03 (curve 
1), 7 = 0.0314 (curve 2) and 7 = 0.05 (curve 3). Panel (c) 
shows the effective potential U(X) for /ii = —0.4, /12 = 0.6 
and different values of a: a = 0.1 (curve 1), a = 1 (curve 2), 
a — 2 (curve 3), a = 3 (curve 4) and a = 5 (curve 5). Panel 
(d) shows the effective potential U(X) for fi\ = 0.2, jj.2 = —0.8 
and different values of a: a = 1 (curve 1), a = 1.7 (curve 2), 
a = 2.5 (curve 3) and a = 4 (curve 4). 



same figure) so that each impurity interacts individually 
with the fluxon. 

Consider now the case /Ui > 0, [12 < 0. The phase por- 
traits for the fluxon dynamics are shown in Figure[7] (panel 
b). The dependence of the threshold pinning current on 
the distance between the impurities is shown in Figure [8] 
(panel b). In the case \x\ = —\i<z = 0.5 at a = the impu- 
rities cancel each other. When a increases, jthr monoton- 
ically increases, tending to the threshold value of one iso- 
lated microshort with the amplitude [i = 0.5. In this case 
trapping occurs only on the microshort because analysis of 
equations (HI)-© shows that 7 t h r (0.5; 1) > "fthr{— 0-5; 1). 
If /ii > |/^2 1 the dependence of j t hr on a is also monotonic. 
At a ~ the fluxon "feels" both impurities as one mi- 
croshort with /zi — \/J-2\- When the impurities separate, the 
contribution of the microresistor to the total amplitude 
weakens and the threshold current gradually increases till 
the value 7t/ ir (Mi> !)■ If Mi decreases and | /X3 1 increases, be- 
haviour of the critical pinning current on a becomes more 
complicated. Consider first the case fi% = 0.4, /i2 = —0.6. 
In the neighbourhood of a = the system can be con- 
sidered as a microresistor with the amplitude \x\ — } . 
When a increases, the well and the barrier, created by the 
microshort /ii start to separate, increasing the depth of 
the well (created by the microresistor). After some value 
of the distance a the dependence of jthr = Jthr(a) ex- 
periences sharp breaking and jthr starts to grow with a. 
The difference between trapping before this breaking point 
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Fig. 8. Threshold pinning current as a function of the distance 
between two impurities ^1/12 < with a = 0.1. Panel (a): 
fj,! = -0.2, fi 2 = 0.8 (A); fix = -0.4, fx 2 = 0.6 (x); fn = -0.5, 
= 0.5 (o); Ml = -0.6, fi 2 = 0.4 (+); a x = -0.8, ju 2 = 
0.2 (V) and m = -0.9, fi 2 = 0.1 (o). Panel (b): jm = 0.8, 
^2 = -0.2 (A); /xi = 0.5, a 2 = -0.5 (o); m = 0.4, = -0.6 
(+); mi = 0.3, a 2 = -0.7 (x); m = 0.2, M2 = -0.8 (V) 
and /ii = 0.1, fi2 = —0.9 (o). Dashed lines correspond to the 
values 7th r (±l; 1). Inset shows more detailed picture of "/thr = 

and after it is based on the trapping scenario at 7 < jthr ■ 
For values of a below the breaking point trapping occurs 
in the well, created by the microresistor (curve 2 of Fig- 
ure [TJd) while after the breaking point trapping occurs on 
the microshort. In other words, the breaking point sig- 
nals the value of the separation of the impurities, before 
which the fluxon "feels" them as one microresistor and 
after which the fluxon "feels" them separately. Decrease 
of fi\, and subsequent increase of I//2I leads to the grad- 
ual shift of the breaking point to the right and smoothing 
of the shape of the dependence jthr{a)- Further decrease 
of (Xi and increase of | /X2 1 makes the dependence 7t/i r (o) 
more and more flat, so that in the limit H\ — ► 1, fi2 — v 
it tends to the horizontal line "fthr — Jthr(~ U !)■ 

4 Conclusions 

We have investigated the fluxon transmission in a dc- 
biased long Josephson junction (LJJ) through two or more 
impurities of different polarity: microshorts and/or mi- 
croresistors. We have observed that the threshold pinning 
current can depend on the distance between impurities in 
the resonant way for the case of two or more microshorts 
or two or more microresistors. That means that at some 
value of the distance the threshold current attains a min- 
imal value, which is less than the threshold current of the 
strongest impurity. The resonant transmission does not 
occur if the fluxon interacts with two impurities of differ- 
ent sign: a microshort and a microresistor. 

The observed effect should not be confused with the 
resonant soliton transmission in the non-dissipative cases 



5,6. In the case of fluxon dynamics in a long Josephson 
junction the presence of dissipation is unavoidable. Far 
away from the impurities fluxon exists as an only one at- 
tractor of the system with the velocity, predefined by the 
damping parameter and external bias. Therefore, contrary 
to the non-dissipative case, there is no sense in comput- 
ing the transmission ratio, which in our case can take only 
two values: zero (trapping) and unity (transmission) . Also 
it should not be confused with the fluxon tunneling as a 
quantum-mechanical object [11] across the double-barrier 
potential, created by two identical microshorts. 

The discussed phenomenon can be observed experi- 
mentally in an annular LJJ via monitoring the current- 
voltage (IV) characteristics. For a LJJ with one impurity 
the fluxon IV curve has a hysteresis-like nature with two 
critical values of the dc bias (discussed in Section^. The 
lower one is the threshold pinning current, which is the 
smallest current for which fluxon can propagate. Although 
simulation in this paper have been performed for the in- 
finite junction, there should be no principal differences 
with the case of an annular junction with sufficiently large 
length L Xj. Currently experiments are performed in 
the junctions L ~ 10Aj (see [2]) that can be considered 
as long. 

For the future research in this direction it is of interest 
to find out how the resonant fluxon transmission changes 
if the actual size of the impurities and the junction width 
along the y-axis are taken into account. 
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